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> ■ ABSTRACT 

OO \ The Parkes Multibeam Pulsar Survey is the most successful survey of the Galactic 

04 ■ plane ever performed, finding over 600 pulsars in the initial processing. We report on 

reprocessing of all 40,000 beams with a number of algorithms, including conventional 
frequency-domain searches and an acceleration search for fast binary pulsars. The 
very large volume of results coupled with the need to distinguish new candidates 
from known pulsars and their many harmonics, often with multiple detections from 
different search algorithms, necessitated the development of a new graphical selection 

■ tool tightly linked to a web-based results database. We discuss and demonstrate the 

benefits of these software systems which are specifically designed for large survey 

Q i projects. The results of this processing have been encouraging; we have discovered 

$-H ' 128 new pulsars including 11 binary and 15 millisecond pulsars, in addition to those 

\ previously found in the survey, we have thus far discovered 737 pulsars. In this paper 

d • we discuss the discoveries of PSR J1744— 3922, a 172 ms mildly recycled pulsar in a 

^ ' 4.6 hr orbit that exhibits nulling behaviour, not previously observed in recycled or 

• i-H , binary objects; PSR J1802— 2124, an intermediate mass binary pulsar (IMBP) and 

' PSR J1801 — 1417, a solitary millisecond pulsar. 

H ' 

' Key words: pulsars: general — pulsars: searches — pulsars: timing — pulsars: indi- 

vidual PSR J1744-3922 — pulsars: individual PSR J1801-1417 — pulsars: individual 
PSR J1802-2124 
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1 INTRODUCTION 

Binary pulsars provide the only environment for precise 
measurements to test general relativity in strong-field condi- 
tions (Taylor 1993). The orbital decay of double neutron star 
systems is currently the only observational demonstration 
of the existence of gravitational waves (Taylor & Weisberg 
1989). For the elTects of general relativity to be seen, the 
pulsar usually has to have a massive companion (> 1 M0) 
and an orbital period of only a few hours. The discovery 
of a black hole - pulsar binary would be very important 
for even more stringent tests of general relativity (Damour 
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& Esposito-Farese 1998). The pulsar itself must be timed 
to high precision, implying that ideally it will be a strong 
millisecond pulsar with a narrow pulse profile. 

Millisecond pulsars with spin periods less than 2 ms 
probe the equation of state for the degenerate matter which 
makes up neutron stars. PSR B 1937-1-21 has the shortest- 
known period of 1.55 ms (Backer et al. 1982). Finding a 
new pulsar with a shorter spin period would clearly be very 
significant. Conversely, not finding such a system in a ma- 
jor survey reinforces the so-called stiff theories of equation 
of state (e.g. Friedman 1995) or other spin-period limita- 
tions, such as gravitational radiation (Andersson, Kokkotas 
& Schutz 1999). 

Finding binary pulsars with orbital periods of a few 
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hours or less is very difBcult; the received signal from pulsars 
in binary systems has a period variation dependent upon 
the orbital motion, making conventional spectral analysis 
considerably less sensitive. To recover sensitivity to binary 
pulsars, particularly in relatively long observations, requires 
the use of computationally expensive 'acceleration' searches 
(e.g. Camilo et al. 2000a). 

We report here the successful use of an efficient accel- 
eration search algorithm and optimised interference filters 
applied to the large number of observations taken as part 
of the Parkes Multibeam Pulsar Survey (PMPS), described 
in Manchester et al. (2001), hereafter 'Paper I'. Section 2 
briefly summarizes the PMPS, with a review of all employed 
search algorithms in Section 3. Section 4 details the data 
reduction processing and Section 5 describes a graphical ap- 
proach to reviewing the results and ensuring that known pul- 
sars were correctly identified. In Section 6 we demonstrate 
that the system performed as oxpoctod. Also, wo report the 
discovery of a significant number of millisecond pulsars us- 
ing conventional processing techniques. These discoveries are 
discussed in Section 7 and future activities are considered in 
Section 8. 



2 PARKES MULTIBEAM PULSAR SURVEY 
SYSTEM 

The PMPS covers the Galactic plane in the region |6| < 5° 
and 260° < I < 50°. It has been very successful, with re- 
ported discoveries of more than 600 new pulsars so far, dis- 
cussed in Paper I, Morris et al. (2002), Kramer et al. (2003) 
and Hobbs et al. (2004). However, there were only four mil- 
lisecond (which we define as having a period < 30 ms) and 
twelve binary pulsars reported. In 2002 a new Beowulf com- 
puter cluster at Jodrell Bank Observatory, COBRA, started 
to be used to process all the data from the PMPS, including 
some data that had never been processed before, employ- 
ing newly developed algorithms (hereafter 'COBRA process- 
ing')- 

The details of the receiving system and its performance 
can be found in Paper I. It is the most sensitive survey of the 
Galactic plane ever undertaken, consisting of 35-min integra- 
tions performed with a sensitive, wide bandwidth receiver. 
Most of the observations took place between 1997 and 2002. 
The data now represent a unique archive, which may not 
easily be reproduced due to steadily increasing radio fre- 
quency interference. By processing the data in a consistent 
fashion, all pulsars above the limiting fiux density thresholds 
within the survey area can readily be compared. The raw (as 
recorded at the telescope) data, the processing results and 
plots are archived using on-line disk storage. These data arc 
expected to provide a basis for further analyses (e.g. detailed 
population studies of Galactic pulsars). 

There arc ~14% more observations than the 2670 point- 
ings of 13 beams originally plarmed. Additional observations 
were taken for a number of reasons including aborted obser- 
vations due to high winds, difficulty in reading some tapes 
thus losing some data, or administrative errors. The ad- 
ditional observations can bo useful duplicates of the main 
observations; for example they can be used for initial val- 
idation of a candidate or a discovery of a pulsar that has 
an intermittent received signal, possibly because of nulling 



(see Backer et al., 1970) or scintillation (see Lyne & Rick- 
ett, 1968). All these data, including truncated observations 
(17.5 - 35 mins), have now been processed. 



3 SEARCH ALGORITHMS 

The standard frequency domain search using a fast Fourier 
transform (FFT) is extremely good at detecting stable peri- 
odicities from solitary pulsars. The standard FFT approach 
was applied to all the data of the PMPS as described in Pa- 
per I. Harmonic summing (Taylor & Huguenin 1969) was 
used to compensate for the relatively low duty cycle of 
a pulsar to maximise sensitivity. However, there are some 
classes of pulsars where the FFT-based search is less sensi- 
tive, including: pulsars in binary systems, 'long-period' pul- 
sars which we define as having periods > 3 s and pulsars 
in which we can observe only intermittent individual pulses. 
Throughout this paper we use the term 'detection' for any 
signals identified by a search algorithm, including known 
pulsars or possible new pulsars and their harmonics and in- 
terference. The term 'candidate' is a detection that has been 
reviewed and considered to be a possible new pulsar. 

3.1 Acceleration Searches 

Searching for pulsars in binary systems is made difficult be- 
cause their changing line-of-sight velocity causes a varying 
received pulsar period. The effect on a search is to spread 
the signal detection over a frequency range. An FFT divides 
the received signal into discrete frequencies, so-called 'spec- 
tral bins', hence, a varying received frequency will spread 
a detection over a rmmbcr of bins, thus reducing signifi- 
cantly the signal-to- noise ratio (S/N) of a detection. This 
effect becomes more severe with shorter orbital and pulsar 
periods, and with longer observations. Consequently, some 
of the most interesting systems, viz, millisecond pulsars in 
tight binary orbits, are very difficult to detect. 

It is possible to modify the time series such that the 
pulse period variations caused by binary motion are re- 
moved; an FFT-based search can then be used to find pe- 
riodicities with no loss of sensitivity. This is a 'coherent' 
search, meaning that all the available information (time, 
amplitude and phase) from an observation is preserved and 
operated on as a whole. However, this method is extremely 
computationally expensive in a blind search, with three pa- 
rameters (binary period, orbital amplitude and phase) to 
search for circular orbits, with an additional two (orbital 
eccentricity and longitude of periastron) for significantly el- 
liptical orbits. This is impractical using today's technology 
(e.g. Dhurander & Vecchio 2001). 

The task is simplified by only searching over one param- 
eter, average acceleration, which approximates the period 
changes over time to a straight line, referred to as a 'lin- 
ear acceleration' search. This works well for orbital periods 
which are much longer than the observation length, partic- 
ularly if the observations happen to be made at an orbital 
phase where the pulsar's acceleration varies least from being 
constant (Johnston & Kulkarni 1991). Linear acceleration 
searches can be performed either in the time or frequency 
domain. In the time domain, the time series is adjusted by 
resampling to compensate for the quadratic change in pulse 
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Figure 1. An illustration of the stack search approach to ac- 
celeration searching. The received time series is first broken into 
a number of segments (in this case 16), which are individually 
FFT'd and the amplitude spectrum formed, using interpolation 
to recover signals lying midway between bins. These spectra are 
then summed at different rates of change of frequency to produce 
a final spectrum. Note that higher harmonics of signals need to 
have a rate of change of frequency proportional to the harmonic 
number. 



phase over the observation time, and then an FFT based 
search is performed. This needs to be repeated for all the 
different accelerations tested as discussed in Camilo et al. 
(2000a). 

Alternatively, in the frequency domain, an FFT of the 
zero acceleration time series is performed and then a series 
of matched filters is applied to correct the frequency spec- 
trum to look for accelerated periodicities (e.g. Ransom et 
al. 2002). While the frequency domain approach is compu- 
tationally more efficient for a blind search, both techniques 
still require considerable processing at each dispersion mea- 
sure (DM) to be searched. Consequently, a coherent linear 
acceleration search is generally only feasible if the DM is 
already known, as in the case of globular cluster searches 
when there is already a confirmed pulsar, ft is currently im- 
practical to apply a coherent acceleration search using 200 
trials (see Camilo et al. 2000a) to the PMPS data, since the 
total processing time would increase by ^100 times. 

By operating on the data in sections, known as an 'inco- 



herent' search, some information is lost, in this case phase. 
However, by sacrificing some sensitivity, an incoherent accel- 
eration search, or 'stack search' as discussed by Wood et al. 
(1991) can be used to save substantial computing resources. 
The stack search is illustrated in Figure 1. It involves split- 
ting the received time series into a number of segments (16 
in the case of the PMPS), FFT'ing each of these individu- 
ally and producing an amplitude spectrum. Use of very short 
segment observations, over time Tseg, considerably improves 
the relative sensitivity to varying periods. In a short obser- 
vation the period does not change much and each spectral 
bin covers a much larger frequency range (A/ = l/Tsog). 
By summing each of the bins of the same frequency from 
all the segments, a spectrum of the full observation can be 
made. Some overall sensitivity will be lost due to the inco- 
herent summing. However, the pulse duty cycle, hence the 
summing of harmonics makes it hard to predict the precise 
effect on spectral S/N. In the case of the PMPS, there are 5 
trials, the fundamental only and then progressively adding 
up to the 2nd, 4th, 8th and 16th harmonics. By compar- 
ing the spectral S/N for the standard search with spectral 
S /N for the stack search for a variety of solitary pulsars we 
found there to be a S/N loss of between 15-35% over a co- 
herent search. Linear accelerations, a, can then be applied 
to the whole sequence by adding up the bins at a slope pro- 
portional to the change in frequency over the observation. 
The frequency change, Av, from segment to segment for an 
observation of length T with number of segments A'sog and 
speed of hght c is: 



Au = voaT / {Nsogc) Hz. 



(1) 



Since the number of bins shifted must be an integer, the fre- 
quency change needs to be rounded to the nearest bin. Prior 
to harmonic summing, the nth harmonic must be corrected 
for acceleration using a slope n times the fundamental slope. 
The FFTs for all of the segments take a similar or less time 
to produce than a single full length FFT. Spectra for each of 
the trial accelerations are produced by addition, making the 
stack search fast and efficient. This is the technique that has 
been applied to processing of the PMPS. Without the stack 
search, the double neutron star PSR J1756- 22t (Faulkner et 
al. in prep.) would not have been detected. 

Highly accelerated systems which have more than 1.5 
orbits over the length of the observation can be found using 
a 'phase modulation' search (Jouteux et al. 2002, Ransom, 
Cordes & Eikenberry 2003). In the case of the PMPS this 
would apply to orbital periods less than 20 min. The algo- 
rithm finds the orbital periodicity by identifying sidebands 
in the power spectrum of the time-series. Any sidebands 
present are summed by stepping through the power spec- 
trum taking short FFTs. This is a fast process and should 
reveal the pulsar period and orbital period. It has been in- 
corporated into the search algorithms and any results will 
be reported in a separate paper. 



3.2 Other search algorithms incorporated 

Frequency domain searches using FFT's are effective for 
most of the range of pulse periods, from the Nyquist limit 
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up to approximately 3 s in the case of the PMPS. How- 
ever, when dealing with periods which are greater than 3 s, 
which lie in the first few thousand spectral bins, the pe- 
riod resolution becomes very coarse. In addition, the pres- 
ence of low-frcqucncy noise in the Fourier spectrum of the 
radio-astronomical time series can be very difficult to re- 
move. By using a time domain search the period resolution 
is much finer and resolves the true spin period more accu- 
rately. A time domain search works effectively on arbitrary 
pulse profiles, therefore is more effective on the very narrow 
pulses typical of long period pulsars. Since the search is for 
long periods, greater than 3 s, the sampling time can be 
considerably increased with no significant loss of sensitivity, 
thus reducing the processing required. An efBcient algorithm 
(Staelin 1969) has been included in the COBRA processing; 
details of this implementation will be described in a sepa- 
rate paper. This algorithm has been successful at finding a 
new pulsar with a 7.7 s period, which will be published in a 
future paper. 

It has been shown (McLaughlin & Cordes 2003) that 
some pulsars can only be detected by observing their indi- 
vidual giant pulses. It is possible to detect these individual 
pulses if they are large enough. A 'single pulse search' for 
dispersed irregular pulses has been incorporated into the 
search routines. Over 30% of all pulsars (new and known) 
detected in the PMPS were also detected in the single-pulse 
search. In addition, the single-pulse search has detected four 
pulsars which were not detectable in the standard periodic 
processing. The results from this search will be presented in 
a future paper. 



4 PROCESSING 

Details of the observing and search analysis systems are 
mostly as reported in Paper I. Modifications and enhance- 
ments are discussed below. 

All the data were gathered at Parkes and stored on 160 
DLT tapes. The tapes have been collected at Jodrell Bank 
Observatory for processing. With the ongoing reduction in 
disk storage costs, all the raw data have been transferred to 
STB of on-line RAID disk systems. Thus, alternative pro- 
cessing approaches can easily be tried in the future. 

The data flow from observation data through processing 
to produce detection lists and plot files is shown in Figure 
2. The processing task is divided into two principal areas: 
separating the data into individual beams and searching for 
major interference signals, followed by detailed searches of 
individual beams. 

The first parts of the processing, stages one and two 
described in Paper I, take the raw data files, split them into 
thirteen individual beams, and perform interference search 
and excision. The resulting datafiles, header files and ac- 
counting files are passed to the search stage. 

During the period of PMPS observations and on-going 
processing, progressively more frequency filters were added 
to combat interference. However, they also restricted the 
amount of spectrum left for searching, particularly at mil- 
lisecond pulsar frequencies. Significant rationalisation of 
these 'interference' frequency bands weis made prior to 
this COBRA processing, which substantially increased the 
amount of spectrum available to find new pulsars. There is 



Table 1. Frequency filters for removing identified interference, 
probably related to mains 50 Hz plus harmonics and local and 
instrumental sources, v is the central frequency of the individual 
filter, i and j are integers and fnyq is the Nyquist frequency of 
the PMPS. In all cases u wraps around the Nyquist frequency 
and returns to lower froc]uonrios. 



Central frequency u 


Width 


Ranges 


Limit 


(Hz) 




a i) 


(Hzi 


X (H- 2J-1) 
8.r92 ^ 


0.001 


\j:l—>W 


V < I^nyq 




0.1 
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\ui xi + U2X j\ 


0.004 


ri:1^50 
tj:-20^20 


V > 1.00 


i/i X i 


0.04 


51 160 




ui X i 


0.4 


1 — > i^nyq4/l'l 




uixi + 3.9 


0.3 


1 Vnyqi/vi 




uixi- 3.9 


0.3 


1 V^yofll^X 




50 X i 


0.2 


1 — > t'nyq 




200 + i 


0.26 


-20 -» 20 




0.14704 X i 


0.0024 


1 ->■ 100 




1349 


12 






508.65 X i 


0.400 X i 


1 -+ 2 




2.045 


0.002 






1.2032 X i 


0.001 X i 


1 ^ 40 




1.0000 X i 


0.001 X i 


1 ^ 30 




0.6433 X i 


0.001 X i 


1 -+ 50 




Frequency definitions: 


VI = 134.99776 Hz, U2 


= 99.995921 


i/nyq = 2000 Hz, 


i^nyq2 — 


4000 Hz, i/nyq4 = 8000 H 



now < 10% of the spectrum removed, as opposed to ~40% 
in previous analyses (Hobbs 2002). This is a significant con- 
tributor to the success at finding a high proportion of new 
millisecond pulsars. A summary of the filters is shown in 
Table 1; a full discussion of the purpose and changes made 
to each of the filters is given by Hobbs (2002). 

Terrestrial interference is normally not dispersed, so an 
FFT is taken of the time series not adjusted for dispersion, 
the 'zcro-DM' spectrum, and the result is stored. These data 
can be used for a future analysis of interference signals at 
Parkes. 



4.1 Beam processing 

The searching of each of the 40,000 beams (stage 3 of Pa- 
per I) can be done on a beam-by-beam basis, essentially 
unafi^ected by any other data or process. This is the perfect 
environment for a Beowulf cluster, which in its simplest form 
is a collection of interconnected PC's. Hence, each beam is 
assigned to a separate processor; each node consisting of 
two processors is able to search two beams simultaneously. 
An individual beam was typically processed in about 6—7 
hours on a single processor. Excellent total throughput of 
processing is achieved by assigning many processors, typi- 
cally about 100. A tape of ~300 beams was processed in a 
day. 

The various tasks (see Figure 2) that make up a search 
of a beam are controlled by a script running on each proces- 
sor. 

The data are de-dispersed as described in Paper I. The 
search routines then process the time-series for each of the 
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325 different DMs. Tfie spectra from each DM for both accel- 
erated and standard FFT searches are searched for the top 
100 detections (increased from 50 in Paper I), which ensures 
that all detections above the noise floor are stored. The re- 
sulting detection lists are stored. To qualify the detections, 
they are analysed for the highest S /N signals (reconstructed 
by inverse transformation to the time domain for standard 
searches and spectral S/N for the segmented search) and the 
DM of the strongest detection is found for each identified 
period. Detections with harmonic relationships are identi- 
fied and eliminated except the period corresponding to the 
strongest signal and the detection with twice this period. 
Any detections over a S/N of 6 are then further analysed in 
the time domain largely as described in Paper I. The stan- 
dard FFT detections are searched over DM and a narrow 
period range, using 64 sub-integrations (increased from 16 
sub-integrations in Paper I) each of ~30 s, to find an opti- 
mised S/N and the results recorded, known as a 'candidate 
plot'. The S/N is calculated by convolving the profile with 
a box-car function, trialed with widths from 1 bin to half 
the length of the profile. Similarly, the best detections from 
the accelerated search are searched in period and frequency 
derivative by repeated adjustment of the time series, and the 
results also recorded. Examples of both types of candidate 
plot are shown in Figure 3. 

For the long-period search, six results from each beam 
which are the best S/N detections from any DM searched, 
are recorded. These are post-processed into candidate plots 
and stored for later inspection. The single pulse search stores 
S/N, time index in the observation, width and DM for all sin- 
gle pulses with a S/N in excess of 4. The period, orbital pe- 
riod and DM information from the phase modulation search 
are recorded for off-line analysis. 



5 FINDING CANDIDATES 

The processing produced a large quantity of output to be 
analysed, with each tape of twenty six pointings producing 
around 50,000 qualified detections. These detections are a 
mixture of interference signals, known pulsars with their in- 
teger, non-integer harmonics and sub-harmonics and, occa- 
sionally, candidates but also their harmonically related de- 
tections. Many of these results will be duplicated in more 
than one of the standard, stack or long period search algo- 
rithms. 



5.1 Benefits of a graphical interface 

With tens of thousands of signals to choose from, it is clearly 
not feasible to manually inspect all the output plots from 
the search algorithms. Two different philosophies to select 
a more reasonable number of candidates for inspection are: 
(a) an algorithm that selects mainly on S/N; (b) give the ob- 
server enough information in a comprehensible form to en- 
able a selection to be made. The algorithmic approach was 
used originally on the PMPS as discussed in Paper I. Using 
a graphical plot to visually select suitable candidates means 
that interference is more easily avoided and fewer plots need 
to be inspected. Candidates that would not be identified by a 
simple algorithm can be seen graphically. For example, if the 
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Figure 2. Processing on COBRA. There are 3 phases: firstly, 
data from the tape are loaded onto archive RAID; the data are 
then copied to COBRA's local working disk and split into beams 
followed by a search for interference specific to a tape; the beam 
data, headers and accounting files are put into a buffer; the beam 
processing searches an individual beam on one processor, followed 
by candidate identification. 



detection does not have a high S/N, but is clear of interfer- 
ence and the noise floor, which can relatively easily be seen 
on a graphical display, then it would be selected for inspec- 
tion, as with the discovery of PSR J1811-24 shown in Figure 
4. These candidates may be the most interesting discover- 
ies. Other benefits of a graphical approach are the relative 
ease of investigations of possible candidates/known pulsars, 
e.g., detections in adjacent beams, harmonic relationships, 
positions of known pulsars relative to beam centres etc. The 
graphical approach has been used successfully for the Swin- 
burne intermediate-latitude pulsar survey (see Edwards et 
al. 2001). 

5.2 Reaper 

A suite of programs, known as Reaper, for graphically in- 
vestigating the results of the PMPS has been developed. 
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File; PM0168_1iai1 RAJ: 19:15:28.0 DecJ: +15:06:27. Gl: 49.968 Gb: 2.122 Date: 020911 
Centre freq. (Hz): 16.93510975 Centre period (ms); 59.04892349 Centre DM: 171.47 
File start (biks): 1 Spectral s/n: 19-7 Recon s/n: 12.7 BIk length (s) 0.76800 L 
Tsomp (ms); 1.0000 Frchi; 1516.5000 DM factor: 1.0 Full Seq: K3436 PSR J1915+ 

Ref MJD: 52528.36088 BC Ref MJD: 52528.36407 



File: PM0168_1 1811 RAJ: 19:15:28.0 DecJ: +1 6:06:27. Gl: 49.968 Gb: 2.122 Dote: 020911 
Centre freq. (Hz): 16.93534826 Centre period (ms): 59,04809189 Centre DM: 168.25 
File start (bIks): 1 Spectral s/n: 22.0 Centre fdot: 1910.00000 BIk length (s) 0.76800 L 
Tsamp (ms): 1.0000 Frchi: 1516.5000 DM factor: 1.0 Full Seq: K3439 PSR J1915-t 

Ref MJD; 52528.36088 BG Ref MJD: 52528.36407 



-10 10 • 



100 200 300 400 



-0.01 -5x10 ' 



1 

Best prd (ms): 59.04714562 

BC prd (ms): 59.04334447 Err: 0.00006317 

Best frq (Hz): 1 6.935620 

BC frq (Hz): 16.936710 Err: 0.000018 

Best DM: 168.55 Err: 1.61 

Best Width: 16 Best SN: 14.5 



1 

Best prd (ms): 59.04902293 

BC prd (ms): 59.04522166 Err: 0.00002868 

Best frq (Hz): 16.935081 

BC frq (Hz): 16.936171 Err: 0.000008 
Best fdot: 1148.49 Err: 865.36 

Best Width: 8 Best SN: 33.1 





n Reaper - AJF 



5-^u.g-2004 2D:26 



-Aug-2004 20:26 



Figure 3. Example candidate plots compajring standard (left) and acceleration (right) searches, these show the re-detection of PSR 
B1913+16, the first binary pulsar. It can be seen that the acceleration search considerably improves the final S/N and pulse profile.The 
acceleration information presented shows beam and detection information at the top; the top plot shows the S/N variation with small 
changes in the frequency due to acceleration (fdot) and period; centre left is a signal strength plot over pulse period (x-axis) and time 
through the observation (y-axis) of 64 sub-integrations; centre right shows four frequency bands plotted against pulse phase; finally the 
corrected pulse profile is along the bottom. The standard search information is similar to the accelerated plot except the top left plot 
shows S/N variation with period changes and DM, the top right plot shows S/N variation with widely changing DM. 



AH the results from one tape of observations are loaded at 
a time. The flow of data from lists of detections and their 
candidate plots, existing pulsars, previous candidates and 
previously viewed candidates through to observation lists at 
Parkes and information on a web site is shown in Figure 5. 
An illustration screen is shown in Figure 6. Its user interface 
consists of an X-Y plot of a pair of chosen parameters, typ- 
ically period and S/N, with a menu of selectable commands 
on the right. The parameters that may be plotted axe: 

• spin period 

• spectral S/N 

• reconstructed S/N (if available) 

• candidate plot S/N 

• dispersion measure 

• pulse width 



• number of harmonics summed 

• spin frequency 

• spin frequency derivative (if available) 

• Galactic latitude and longitude 

Other information can be superimposed with the detec- 
tions, including known pulsars and their harmonics, previ- 
ously selected candidates and previously viewed detections. 
This is done using a variety of colours and symbols. 

The X-Y plot can be navigated with zoom and move- 
ment commands. Any detection can have its candidate plot 
viewed by selecting its symbol. The Galactic position (see 
Figure 7) of a possible candidate can be displayed. Also on 
the plot, the exact position, period and DIvl of known pul- 
sars can be shown, primarily to confirm that there is no 
potential link between the detection and a known object. 
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Figure 4. Annotated Reaper screen-shot of the discovery of PSR 
J1811-24 showing that the S/N is only slightly above the noise 
floor and surrounded by interference, but still easily discernable. 
Other isolated candidates were seen to be interference from their 
candidate plots. 

Visual guides for the beam width and beam identifications 
can be overlaid, hence showing the beams which may have 
other detections. 

If necessary, all the detections for a beam can be viewed, 
for example, to investigate an apparent non-detection of a 
known pulsar or see a search's response to a very large inter- 
ference signal or powerful nearby pulsar. All the detections 
from the standard and accelerated searches are available. In 
this case the top 100 detections at each DM and harmonic 
sum are recorded with S/N and, if appropriate, the accel- 
eration in terms of frequency derivative. This amounts to 
more than 300,000 detections per beam. Graphically these 
can be viewed and zoomed into for particular periods over 
the various DMs seen. 

The graphical access to the results database for the en- 
tire PMPS is a valuable tool for investigating known pulsars, 
interference signals and potential single-pulse candidates. 



5.3 Selecting and confirming candidates 

The selection of candidates is a matter of experience; in gen- 
eral candidate plots are viewed if they have S/N >8, are not 
identified as the fundamental or harmonic of a known pulsar 
and are clear of detections in multiple beams, which would 
probably be interference. Solitary candidates normally are 
present in both the standard and acceleration searches. If 
the plot shows a reasonably good profile at a DM which is 
greater than 3 pc cm~^ (to avoid terrestrial interference) and 
does not get changed substantially from the spectral search 
through time-domain refinement (see Section 4.1) then it 
may be a candidate. It is important to verify that it is not 
related to a known pulsar, so, a check is made of the im- 
mediate area of sky ±1° using Reaper's Galactic plot. A 
final check is made by restricting a S /N - period plot to just 




Figure 5. Data flow after processing, through Reaper , to cre- 
ate now candidates for observation at Parkes and placing on a 
website. See text for 'Best Obs.' description. 
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Figure 6. A typical search screen for Reaper , showing detec- 
tions, candidates and known pulsars for a particular tape of obser- 
vations. The clickable menu of commands on the right highlights 
chosen selections; the details are not relevant here. 



the detection beam and ensuring that there are no large in- 
terference signals in the beam which could cause spurious 
detections. 

After carefully considering a detection to ensure that 
it is unlikely to be interference or a known pulsar, it may 
be assigned to one of three classes as a candidate. Class 1 
candidates will be prepared for confirmation observations at 
Parkes at the next opportunity; class 2 candidates are kept 
for further consideration and class 3 candidates are effec- 
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Figure 7. Reaper displaying a Galactic plot indicating FWHM 
of the observation beams and highlighting known pulsars. 

tively discarded but kept for future reference. An individual 
detection can be promoted or demoted by reselecting their 
class within Reaper. 

At any time, candidates selected are in various stages of 
the confirmation process - waiting for confirmation, obser- 
vations made but not seen, gridded (e.g Morris et al. 2002), 
confirmed as a new pulsar or rejected. To keep all the in- 
formation readily available. Reaper writes out candidates 
to a web site; other information is also collected from the 
central archive files and displayed. The candidate plot is an 
essential feature, so, confirmation plots are also displayed 
for ease of reviewing at the participating institutions. 

Reaper has been used to view all the processing results 
from the PMPS, typically only taking about one hour per 
tape. Both the processing system and Reaper are modular 
and can be ported to other major surveys. For example, both 
systems have been ported to the related survey, the Parkes 
High Latitude Pulsars Survey (the PH Survey) . It was in the 
PH Survey that Reaper first detected the millisecond pul- 
sar PSR J0737— 3039A of the double-pulsar system (Burgay 
2004). 



6 SYSTEM PERFORMANCE 

With a complex system such as the PMPS, including the 
subsequent processing, it is important to know that it per- 
forms as predicted. To verify if the sensitivity, sky coverage 
and detection algorithms are working correctly, we consider 
the re-detections of known pulsars, pulsars in the search area 
which were not detected and any detections that showed any 
anomalies. 

There are more than 1000 pulsars in the survey area 
found from many surveys including all the PMPS process- 
ing. As has been reported in Hobbs et al. (2004), of the 
264 pulsars in the survey area discovered by other surveys, 
all but eleven have been re-detected with no unaccountable 
non-detections. 

Pulsars were detected over the whole area of the survey 
at spin periods from 1.85 ms to 7.7 s and at DM's from 3 to 
1200 pc cm~^. Reviewing flux density measurements from 



Table 2. Provisional names and parameters for all the binary and 
millisecond pulsars found using COBRA processing, showing spin 
period, orbital period and DM. The spectral S/N for the standard 
search (Std.) and the stack acceleration (Acc.) search are shown 
to highlight the differences in search algorithm. The 'Plot' S/N is 
the refined S/N from the candidate plot, the values in italics are 
where the acceleration plot S/N has been used. 



PSR Jt 


Period 


Orbit 


DM 




S/N 






(ms) 


(days) (pccm ■') 


Std. 


Acc. 


Plot 


1125—60 


2.630 


8.75 


53.1 


11.5 




9.8 


1215-64 


3.540 


4.08 


17.7 


19.5 


16.1 


16.6 


1439-54 


28.639 


2.12 


14.5 


7.1 


10.1 


17.5 


1551-49 


6.284 




114.6 


12.2 




10.2 


1723-28 


1.856 


tbd* 


19.9 


16.7 


40.3 


20.8 


1726-29 


27.081 




60.9 


8.9 


7.3 


12.2 


1744-3922 


172.449 


0.19 


145.7 


15.3 


13.8 


26.5 


1756-22 


28.456 


0.32 


121.6 




10.6 


18.4 


1801-1417 


3.625 




57.2 


31.5 


17.5 


23.7 


1802-2124 


12.643 


0.70 


145.7 


8.3 


15.0 


16.1 


1813-26 


4.430 




122.5 


14.7 


9.4 


11.1 


1822-08 


834.856 


290.2 


165.0 


28.7 


17.0 


28.7 


1841+01 


29.773 


10.48 


125.9 


11.1 


7.1 


11.2 


1843-14 


5.471 




114.6 


12.5 


6.2 


11.6 


1853-1-13 


4.092 


115.7 


30.6 


20.5 


10.7 


15.9 


1911+13 


4.626 




31.1 


29.0 


17.7 


38.2 


1910+12 


4.984 


58.32 


38.1 


20.8 


14.6 


18.0 



f Interim pulsar names have two digits of declination. 
* The orbit for PSR J1723-28 has yet to be determined. 

Hobbs et al. (2004) of the weakest detected pulsars shows 
that the system meets the sensitivity calculated in Paper I. 



6.1 Non-standard detections 

The detection of PSR B1635-45, (Johnston et al. 1992 ) was 
at a fundamental period of 529.1 ms instead of the catalogue 

period of 264.55 ms. Further checking showed that the cat- 
alogue period was incorrect and it has now been updated. 

As noted in Section 4, fllters are used to eliminate 
some known interference signals. These filters can also re- 
move pulsars which happen to be at the same period. In 
the case of PSR B1240— 64, which has a spin frequency of 
2.574Hz, its detection was discarded, with the first 12 har- 
monics, by a filter which removed frequencies from 0. 6423Hz 
to 0. 6443Hz together with the next 50 harmonics. Because 
PSR B1240-64 is very strong, having a S/N of 1466, it was 
seen at its 13th and other higher harmonics and strongly 
at a period of 3/4 of the fundamental; however, a weaker 
pulsar might have been lost. 

Very strong pulsars are seen in multiple beams. One of 
the tests for interference signals is that they are not localised 
to a small section of sky, but seen in many of the beams on a 
tape. These signals are excised from all the beams on a par- 
ticular tape, assuming them to be interference. This did not 
prove to be a problem for the PMPS, since in all these cases 
the pulsar was seen in adjacent beams on a different tape. 
However, at the periphery of a search the very strongest 
pulsars could be lost. 
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6.2 Multiple detections 

Reasonably strong pulsars or those near the edge of beams 
are likely to have detections in adjacent beams. All these 
detections are recorded for analysis and the strongest S/N 
listed as the 'best observation' as shown in Figure 5. 
Reaper can relate beam position to pulsar position and any 
anomalous behaviour can then be identified; for example, 
these detections were used to improve the position of PSR 
J1744— 2335 which is close to the ecliptic plane, where tim- 
ing alone does not give good resolution in ecliptic latitude. 
In this case scintillation was not a problem. Also, nulling 
pulsars can be identified, either by seeing nulling within an 
observation or by finding that a pulsar is not detected in the 
beam(s) in which a detection is expected. 

A particular use of multiple detections is to extend the 
timing baseline of newly discovered pulsars. After a year of 
timing observations it is possible that a phase connected so- 
lution may be found which will be coherent with the original 
detection observation. Including a much earlier pulse arrival 
time will add significant constraints to the cphemeris of the 
new pulsar. Ideally, there will be multiple detections in the 
survey, which can be identified using Reaper to further im- 
prove the timing solution. We note that major surveys us- 
ing multiple beam receivers use interlaced observations to 
have complete coverage of an area. For future timing of new 
discoveries using the survey observations, the scheduling of 
adjacent beams to be evenly spread over the survey period 
is worth considering. 



7 DISCOVERIES 

There have been 128 new pulsars discovered in the COBRA 
processing. These are being regularly observed to obtain a 
coherent timing solution and tlicy all will be published in fu- 
ture papers. A list of provisional parameters for the eleven 
binary and fifteen millisecond (spin period < 30 ms) pulsar 
discoveries axe shown in Table 2. All but two of the binary 
discoveries are also millisecond pulsars. Listed for compar- 
ison arc the spectral S/N for both the standard and stack 
acceleration searches and the S/N from the candidate plot. 
For solitary and long orbital period binaries the S/N for 
the stack search is lower than in the standard search, as ex- 
pected. However, for the short spin period, tight orbit pul- 
sars the S/N with the stack search has been substantially 
improved; as has been stated, for PSR J1756— 22 there was 
no detection at all in the standard search. It is worth noting 
that for very short spin periods (< 10 ms), the spectral S/N 
is typically higher than the candidate plot S/N which was 
derived in the time domain. This is due to the small number 
of time samples over the period of the pulsar, making the 
time domain S/N hard to determine accurately. Reaper was 
used to investigate spectral S/N against period from < 1 ms 
to 20 ms, specifically to look for short-period candidates. 

In Figure 8 we plot spin period, orbital period and DM 
of these new discoveries with the other known millisecond 
and binary pulsars previously found in the search area. It can 
clearly be seen that the COBRA processing, with improved 
interference filters and acceleration searches, has been sub- 
stantially more successful at finding pulsars in binary sys- 
tems with tight orbits and shorter pulsar periods than the 



earlier PMPS processing. These objects are being timed and 
will be discussed in future papers. Paper I discusses the re- 
duced sensitivity of the PMPS to pulsars with very short 
spin periods at increasing DM just due to dispersion smear- 
ing. For example, the sampling time is increased to 1.0 ms 
above a DM of 140 pc cm~^, to match the dispersion smear- 
ing and thus save processing time. It can be seen in Figure 
8(b) that there are no PMPS millisecond pulsars above a 
DM of 150 pc cm"^. 

In the following sub-sections we discuss three newly dis- 
covered pulsars, two in binary systems and a solitary mil- 
lisecond object. Both of the binary systems have very low ec- 
centricity: for such systems, the longitude of periastron is not 
well defined. To avoid large covariance between the usually 
quoted longitude and epoch of periastron, we have adopted 
the ELLl timing model, which uses the Laplace-Lagrange 
parameters rj = esinw and k = ecosw, where e is the or- 
bital eccentricity and lu is the angle of periastron, plus the 
time of ascending node Tasc- A detailed description of this 
binary model can be found in Lange et al. (2001). The pa- 
rameters for the named pulsars are in Table 3, including the 
derived parameters: characteristic age, Tc = P/{2P); surface 
dipolc magnetic field strength, = 3.2 x lO"' (Ppy^^ G. 
The pulsar mass was assumed to be 1.35 Mq. Distance, d, is 
computed from the DM with both the TC93 model (Taylor 
&: Cordes 1993), for consistency with earlier publications, 
and NE2001 (Cordes & Lazio 2002) for comparison. The 
pulse profiles are shown in Figure 9. All of the pulsars re- 
ported here have had regular timing observations covering 
more than one year, either at Parkes or at Jodrell Bank Ob- 
servatory or both. Timing analysis was carried out using the 
TEMPO program^ as described in Paper I. 



7.1 PSR J1744-3922 

PSR J1744-3922 was discovered in February 2003; it com- 
pletes more than 12% of its orbit during a 35 min search 
observation. Consequently, the acceleration correction im- 
proved the optimised S/N by nearly 50% from 17.8 to 26.5. 
There have been regular observations at Parkes since the 
discovery. Subsequently, pulse times of arrival from survey 
observations with detections of PSR J1744— 3922, taken in 
February and June 2001 have been used to improve the 
timing solution substantially. The companion appears to be 
light, with a minimum mass of 0.083 solar masses. This pul- 
sar has an unusual combination of characteristics: (a) it is 
in a compact circular orbit of only 4.6 hours with a pro- 
jected semi-major axis of 0.21 It-sec; (b) it has an unusually 
long spin period of 172 ms and; (c) it experiences a lack of 
detectable emission for some of the time. The spin period 
would suggest a young pulsar, however, the characteristic 
age of 1.8 Gyr and a low eccentricity orbit indicates a mildly 
recycled object. 

The evolution of this type of system is the source of 
considerable discussion in the literature; however a typical 
scenario for a low mass binary pulsar with P), < 1 day is 
given by Phinney & Kulkarni (1994). In summary, the larger, 
primary star with mass greater than the critical mass for a 

t See http://www.atnf.csiro.au/research/pulsar/tempo/ 
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Figure 8. All binary and millisecond pulsars within the PMPS search area. The systems discovered by the PMPS are shown with the 
new discoveries enclosed in a circle. Panel (a) shows the relationship between spin and orbital periods of the binary pulsars. Panel (b) 
shows the spin period and DM for the millisecond pulsars. Named pulsars are discussed in the text. 



supernova explosion (~ 8 Mq) evolves and expands caus- 
ing the secondary to spiral-in. The He core of the primary 
eventually explodes forming a neutron star. During the evo- 
lution of the secondary, mass transfer to the neutron star 
takes place, thus spinning up the pulsar. The second mass 
transfer phase circularises the orbit and probably causes a 
second spiral-in phase, hence, the very tight orbit. 

The position of PSR J1744-3922 is well within the 
95% error box of EGRET source 3EG J1744-3934 (Hart- 
man et al. 1999 ). Indeed, this pulsar was independently 
detected in a targeted search of EGRET sources (Hessels et 
al. 2004). However, since the source is highly variable (vari- 
ability index, V , of 5.12, McLaughlin 2001 ) and sources due 
to pulsars typically have < 1, the association is unlikely. 
Kramer et al. (2003) reviewed possible PMPS discovery as- 
sociations with EGRET 7-ray sources. We follow the same 
assessment approach using the pulsar's spin down energy 
given by: E = Atz^IPP~^ in erg , where a neutron- 
star moment of inertia / = 10*^ g cm^^ is assumed. For 
PSR J1744-3922, E = 19.68 x 10^° erg s"^ at a distance, d, 
of 3.1kpc (using NE2001, Cordes & Lazio 2002) such that 
\og[E/dP (erg s'^ kpc-2)]=30.35. 3EG J1744-3934 has a 
7-ray flux, F, of 1.27 ± 0.26 x 10~^"erg s"^ cm"^ implying 
a 7-ray luminosity which is 550 times the spin-down lumi- 
nosity. Clearly, the EGRET source cannot be identified with 
the spin-powered pulsar PSR J1744-3922. 

The occasional lack of radio emission could be intrin- 
sic or extrinsic to the pulsar. Possible extrinsic phenomena 
are scintillation (e.g. Lyne & Rickett 1968) or the companion 
star obscuring the pulsar. There is no evidence of systematic 
changes in the radiation characteristics at different phases 
of the orbit, so it is unlikely to be due to eclipsing by the 
companion. With a DM of 148 pc cm~"^ and an observing 
frequency of 1.4 GHz, we expect strong difFractive scintillar 
tion to produce scintles with a typical width of a few MHz 



(Cordes et al. 1985). The bandwidth of the receiving system 
is 288 MHz (see Paper I), consequently, any strong scintil- 
lation effects are averaged out and would not be seen. It is 
likely that the cause of the radiation fluctuations is intrinsic 
and probably due to nulling, first reported by Backer (1970) 
for normal pulsars. 

The nulhng of PSR J1744-3922 occurs in short 
timescales. The pulse can be visible from less than a minute 
to a few minutes with off-times up to a few tens of minutes. 
Some 35 min observations do not have any detections. Over- 
all, there is no detectable radiation at 20 cm for ~75% of 
the time. To our knowledge this is the first time that nulling 
has been observed in a binary or a (mildly) recycled pulsar. 

There have been attempts to find correlations between 
pulsar nulling and other pulsar characteristics (e.g. Rankin 
1986 and Biggs 1992); for instance, Biggs found a marginally 
significant correlation between pulsar characteristic age and 
null fraction. For recycled pulsars many parameters are con- 
siderably different e.g. they have much larger characteristic 
ages and their periods are typically shorter. Until this dis- 
covery, recycled pulsars were not observed to null. Physi- 
cally, it is possible that the potential drop, A'l', between 
the magnetic pole and polar cap may have an infiuence on 
a pulsar's emission mechanism. We have investigated the 
magnitude of this potential for different types of pulsars to 
look for evidence of a correlation. In terms of observable 
parameters, the potential drop: 

, -3/2 / p \ 1/2 
A*.2xl0-(-) (^^j V, (2) 

where P is spin period and P the period derivative (Gol- 
dreich & Julian 1969). We considered the median potential 
drops for the groups of pulsars and compared the results 
with PSR J1744-3922. The results are shown in Table 4. 
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Table 3. Positions, flux densities, measured and derived parameters for the newly discovered pulsars with complete timing solutions. 
Values in parenthesis are twice the nominal TEMPO uncertainties in the least significant digits quoted, obtained after scaling TOA 
uncertainties to ensure xi = 1- The pulsar has an assumed mass of I.SSMq. 



Parameter PSR J1744-3922 PSR J1801-1417 PSR J1802-2124 



Right ascension (JiuOU) (n m s) 


17:44:02.675(10) 


18:01:51.0771(3) 


18:02:05.3352(3) 


Ueclmation (JiUUU) ( ) 


-39:22:21.1(4) 


-14:17:34.547(3) 


-21:24:03.6(3) 


Galactic longitude (°) 


350.91 


14.55 


8.38 


Galactic latitude (°) 


-5.15 


-1-4.16 


4-0.61 


Period P (ms) 


172.444360995(2) 


3.6250966597450(15) 


12.647593582763(1 


Period derivative P (x 10~^^) 


1.55(12) 


0.00528(5) 


0.072(1) 


Epoch (MJD) 


52530.0000 


52340.0000 


52855.0000 


Dispersion Measure DM (pc cm~^) 


148.1(7) 


57.2(1) 


149.6(1) 


Orbital Period Pj,(days) 


0.19140624(7) 




0.6988892434(8) 


a sin i (It-s) 


0.2120(2) 




3.718866(7) 


Laplace-Lagrange parameters: 








V 


0.8(20) Xl0~3 




0.3(8) Xl0~5 


K 


0.1(2)xl0~^ 




0.1(6) xlO~^ 


Tasc (MJD) 


52927.16723(3) 




52595.7950781(4) 


Eccentricity e 


<0.006 




<0.00001 


Ntoa 


36 


64 


32 


Timing data span (MJD) 


51953 - 53107 


51561 - 53119 


52605 - 53105 


RMS timing residual (^s) 


197 


26.0 


7.0 


1^ lux density at i4UuMHz (mjyj 


0.20(3] 


0.87(10; 


0.77(9] 


Luminosity, log Sd} (mjy kpc^) 


0.62 


0.45 


0.92 


Wso (ms) 


3.4 


0.35 


0.37 


Wio (ms) 






0.74 


Characteristic Age tc (Gyr) 


1.8 


10.8 


2.8 


Magnetic Field B,, (10*^ G) 


166 


1.4 


9.7 


Minimum companion mass {Mq~) 


0.083 




0.81 


Distance d (kpc) - TG93 


4.6 


1.8 


3.3 


Distance d (kpc) - NE2001 


3.1 


1.5 


2.9 


2 (kpc) 


0.41 


0.13 


0.04 


Detection beam number 


4 


9 


4 


Radial distance (beam radii, ~7 arcmin) 


1.3 


0.8 


0.9 


Detection signal-to- noise ratio (S/N) 


26.5 


23.7 


16.1 



It can be seen that nulling pulsars may have a lower po- 
tential drop than the median of all pulsars or all recycled 
pulsars, which we have somewhat arbitrarily defined as hav- 
ing a characteristic age > lO'^ yrs and surface magnetic field 
< 10" G. There is a difference between recycled pulsars 
with P < 100 ms and those with P > 100 ms, which we 
show separately in Table 4. The lowest potential drop is for 
recycled pulsars with P > 100 ms, which is the category for 
PSR J1744— 3922. With some exceptions, there could be a 
trend among normal pulsars between nulling and A>]/; the 
same relationship between PSR J1744— 3922 and other sim- 
ilar objects is not clear. 



7.2 PSR J1802-2124 

This system has a heavy companion with a minimum mass 
>0.8 Mq which could be another neutron star. However, the 
circular orbit of the pulsar implies that it is unlikely that a 
second supernova explosion occurred in the system, hence 
the companion is probably a heavy white dwarf (e.g. Phin- 
ney & Kulkarni 1994). The evolution is likely to be similar 



Table 4. The potential drop A$ for groups of pulsars. All the 
published data are fallen from the ATNF Catalogue (Manchester 
et al. 2004). The groups used are: [1] all pulsars; [2] pulsars known 
to null taken from Rankin (1986), Biggs (1992) and PMPS pulsars 
observed to null; [3] all recycled pulsars (magnetic field < 10^^ G 
and characteristic age > 10® yrs); [4] short-period recycled pul- 
sars (P < 100 ms); long-period recycled pulsars (P > 100 ms). 
Included is PSR J1744-3922 for comparison. 



Type of Pulsar 


Number 


log[A*(V)] 
Median 




All pulsars 


1381 




[1] 


Known nuUers 


29 




[2] 


All recycled 


121 


14 1+0-6 


[3] 


'short' recycled 


79 


14.3lg:« 


[4] 


'long' recycled 


42 


I3.0l°;| 


[5] 


PSR J1744-3922 


1 


13.0 
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1 1 1 1 

PSR J1744-3922 
0.1724, 148 


1 . J 


. PSR J1801-1417 
A 0.0036, 56 


! 


PSR J1802-2124 
0.0126, 149 

1 1 1 1 



Figure 9. Pulse profiles for the three pulsars discussed. The high- 
est point in the profile is placed at phase 0.3. The pulsar name, 

period (in s) and DM (in pc cm~^) are given. The small horizontal 
bar under the DM indicates the effective resolution of the profile, 
including the effects of interstellar dispersion. Integration times 
used were 313 min, 147 min and 599 min for PSRs J1744-3922, 
J1801-1414 and J1802-2124, respectively. 



to that of PSR J1744— 3922, except that the progenitor sec- 
ondary star was considerably more massive in this system. 
This would result in the neutron star being engulfed during 
the red giant phase of the companion, accreting material 
and consequently being spun-up. This process would also 
cause the two stars to spiral-in towards each other and to 
eject the giant's envelope, leaving a pulsar in a tight, circu- 
lar orbit with a large white dwarf companion. The system 
is categorised as being an intermediate mass binary pulsar 
(IMBP) as introduced by Camilo et al. (1996) and further 
discussed in Camilo et al. (2003), defined as having interme- 
diate mass donor staxs. They identified binary systems with 
the pulsar spin period 10< P <200 ms, e < 0.001, compan- 
ion mass >O.4M0 and 0.5<Pb <15 days as likely IMBP's. 
PSR J1802— 2124 meets all these criteria and is clearly a 
member of this class of system; it has added to the sub-class 
of relatively short spin periods with a 12.6 ms rotation, sim- 
ilar to PSR's J1435-6100/J1757-5322 which contrast with 
PSR B0655+54 at 195 ms noted by Camilo ct al. 

PSR J1802-2124 brings the number of IMBPs known 
to 15 (for a list of the other 14 see Hobbs et al. 2004); with 
eight either found or detected by the PMPS. This concentra- 



tion near to the Galactic plane is discussed in Camilo ct al. 
(2001), with a possible reason being the birth velocity is re- 
stricted because of the high system mass. PSR J1802— 2124 
with \z\ of only 0.04 kpc supports this view, although there 
may bo selection effects due to the high sensitivity survey 
coverage of the PMPS along the Galactic plane. 

PSR J1802-2124 has a narrow pulse profile of < 0.4 ms 
(< 3% of the period), which may be even narrower if the 
effects of dispersion smearing are removed. This probably 
accounts for the flux density being apparently less than the 
nominal PMPS flux limit, which is calculated at a 5% duty 
cycle. 

7.3 PSR J1801-1417 

PSR J1801-1417 is a solitary 3.62 ms pulsar flrst discov- 
ered in November 2002. It has a characteristic ago of 10.8 
Gyr. The ephemeris of this pulsar has been considerably im- 
proved by including three detections from the PMPS back 
to January 2000. 

Following an earlier suggestion by Bailes et al. (1997), 
Kramer et al. (1998) put forward evidence that solitary mil- 
lisecond pulsars are less luminous than when they are in 
binary systems. By calculating the flux density times the 
square of the distance, ScP, we can make luminosity com- 
parisons; the values of log Scf (rn,]y kpc^) for our sources 
are shown in Table 3. As can be seen the luminosity of 
PSR J1801-1417 is significantly lower than that of PSR 
J1802— 2124, although both are more luminous than the 
ranges identified by Kramer et al. (mean value, log Sci^ of 
—0.5 ± 0.3 and 0.2 ± 0.1 for solitary and binary millisecond 
pulsars respectively). We note that it is unlikely that either 
of these pulsars would have been discovered if their lumi- 
nosity was as low as Kramer et al.'s ranges, both pulsars 
exceed the 1.5 kpc previously used, with PSR J1802— 2124 
at 3.3 kpc. 

We have updated the luminosity statistics for recycled 
pulsars, as defined by Kramer et al. (1998) for the current 
known population using a distance limit of 2.0 kpc, increased 
from 1.5 kpc. The population has increased from 7 to 10 
solitary pulsars and from 11 to 19 in binary systems. With 
a distance limit of 2.0 kpc the luminosity mean values are 
log Sep of -0.4 ± 0.2 and 0.1 ± 0.1 for solitary and binary 
objects respectively. These results, which were also checked 
at a limit of 1.5 kpc, are consistent with Kramer ct al. and 
provide further evidence that solitary millisecond pulsars are 
less luminous radio sources than their binary counterparts. 
The cause of this effect is still not known. There will be 
future analysis of all the new millisecond discoveries of the 
PMPS. 



8 FUTURE WORK 

The PMPS has been the most successful search for pulsars 
ever undertaken. However, there is still further work to do. 
The processing performed so far is still insensitive to fast 
millisecond pulsars in binary orbits of a few hours or less; 
for example the binary pulsar PSR J0737— 3039 (Burgay et 
al. 2003), would only be detected at the times in its 2.4 hour 
orbit where the apparent spin period changes are most con- 
stant. This system was found using essentially the same soft- 
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ware system as used in this processing, but the observation 
time was only four minutes, which radically reduces pulse 
smearing. Short observations also reduce potential sensitiv- 
ity, so, there is considerable scope for further reanalysis of 
the data with more sophisticated software or as partial ob- 
servations. 

We are processing the observations in four quarters, us- 
ing the standard and stack search algorithms. With a shorter 
observation time there will be less pulse smearing of highly 
accelerated pulsars as discussed in Section 3.1. This may 
find objects in fast orbits, but at a price of reduced base 
sensitivity by a factor of 2. Evolutionary population studies 
suggest that black hole - pulsar binaries exist in the Galactic 
field (Sigurdsson 2003). Such systems may be in very tight 
orbits. One or more could be in the PMPS dataset and a 
better acceleration search algorithm may find them. 
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